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.<^ Abstract: One of the main implications of the LHC discovery of a Higgs boson with a 

mass Mh ~ 126 GeV is that the scale of supersymmetry-breaking in the Minimal Super- 
symmetric Standard Model (MSSM) might be rather high, Ms ^ Mz- In this paper, we 
r~| consider the high Ms regime and study the spectrum of the extended Higgs sector of the 

^^ MSSM, including the LHC constraints on the mass and the rates of the observed light h 

C^ state. In particular, we show that in a simplified model that approximates the important 

r~| radiative corrections, the unknown scale Ms (and some other leading SUSY parameters) 

can be traded against the measured value of M/j. One would be then essentially left with 
only two free parameters to describe the Higgs sector, tan/3 and the pseudoscalar Higgs 
mass Ma, even at higher orders. The main phenomenological consequence of these high 
OO Ms values is to reopen the low tan/3 region, tan/^ <3-5, which was for a long time buried 

under the LEP constraint on the lightest h mass when a low SUSY scale was assumed. 
We show that, in this case, the heavier MSSM neutral H/A and charged H states can 
Ci be searched for in a variety of interesting final states such as decays into gauge and lighter 

^ ; Higgs bosons (in pairs on in mixed states) and decays into heavy top quarks. Examples of 

L| sensitivity on the [tan/3. Ma] parameter space at the LHC in these channels are given. 
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1. Introduction 

The discovery of a Higgs-like particle by the ATLAS and CMS collaborations [1] in July 
2012 was a triumph for the Standard Model (SM) of particle physics as it led to a first 
verification of one of its cornerstones, the Higgs sector [2-4] that realises the breaking of the 
electroweak symmetry and generates the masses of the fundamental particle masses. This 
discovery had also very important consequences on theories beyond the SM, among which 
supersymmetric theories (SUSY) stood as the most promising ones. This is particularly 
the case of their minimal low energy realization, the Minimal Supersymmetric Standard 



Model (MSSM) [5,6] in which the electroweak symmetry breaking sector is extended to 
contain two Higgs doublets which, after symmetry breaking, lead to the existence of five 
physical states: two CP-even h and H, one CP-odd A and two charged H^ bosons [3, 7]. 

The Higgs observation at the LHC with a mass of approximately 126 GeV first gave 
support to the MSSM in which the lightest CP-even h boson was predicted to have a mass 
less than ~ 130 GeV [8]. An annoying problem is that the measured mass value is too close 
to the predicted upper limit on Mh in the MSSM, suggesting that the SUSY scale is rather 
high, Ms > 1 TeV; see for instance the discussion of Ref. [9]. The fact that Ms is large is 
backed up by direct SUSY particle searches, which set limits of the order of 1 TeV for the 
strongly interacting superparticles [10]. In addition, with the precision measurements of its 
couplings to fermions and gauge bosons, the Higgs state looked more and more SM-like, 
as no significant deviations from the SM expectation is presently observed [10]. Although 
this had to be expected since, as is the case in many extended Higgs sectors, there is a 
decoupling limit [11] in which all the heavier Higgs particles decouple from the SM spectrum 
and one is left only with the lightest h state which has almost the SM properties, this is 
again unfortunate. Tests of the properties of the observed Higgs state have to be pursued 
with more accuracy in order to pin down small deviations from the SM prediction. 

An equally important way to probe the MSSM is to search for the direct manifestation 
of the heavier H, A and H states. These searches are presently conducted by the ATLAS 
and CMS collaborations in the regime where tan/3, the ratio of the vacuum expectations 
values of the two Higgs fields, is very large, tan/3 = ^2/^1 ^ 5-10, which significantly 
enhances the Higgs production rates at the LHC. The regime with low tan/3, tan/3 < 3-5, 
is ignored, the main reason being that if the SUSY scale should not exceed Ms ~ 3 TeV 
to have a still acceptable fine-tuning in the model [12], the h mass is too low and does not 
match the observed value. More precisely, this tan/3 region was excluded by the negative 
Higgs searches that were performed at the ancestor of the LHC, the LEP collider [13]. 

In this paper, we reopen this low tan /3 region by simply relaxing the usual assumption 
that the SUSY scale should be in the vicinity of 1 TeV. In fact, many scenarios with a 
very large scale Ms have been considered in the recent years, the most popular ones being 
split-SUSY [14] and high-scale SUSY [15]. In these constructions, the SUSY solution to 
the hierarchy problem is abandoned and the masses of all the scalars of the theory (and 
eventually also those of the spin-^ superparticles in high-scale SUSY) are set to very high 
values, Ms ^ Mz- Hence, all the sfermions and Higgs bosons are very heavy, except for 
a light SM-like Higgs boson whose mass can be as low as M^ ~ 120 GeV even if tan/3 is 
very close to unity. In fact, for this purpose, the scale Ms needs not to be extremely high, 
for instance close to the unification scale as in the original scenarii of Refs. [14, 15], and 
values of Ms of order 10 to 100 TeV would be sufficient. 

In addition, one may assume that only the sfermions are very heavy and not the Higgs 
particles, as it is would be the case in non-universal Higgs models where the soft-SUSY 
breaking mass parameters for the sfermion and the two Higgs doublet fields are disconnected 
[16]. On would have then a scenario in which the entire MSSM Higgs sector is kept at the 
electroweak scale, while the sfermions are pushed to the high scale. Such scenarios are also 
being considered [17] and they might occur in many theoretical constructions. 



A first important aspect that we will address in this paper is the treatment of the 
radiative corrections in the Higgs sector and the derivation of the superparticle and Higgs 
spectrum in these high scale scenarios. It is well known that for Ms values in the multi- 
TeV range, the MSSM spectrum cannot be obtained in a reliable way using the usual RGE 
codes that incorporate the higher order effects [18, 19]: one has first to decouple properly 
the heavy particles and to resum the large logarithmic contributions. Such a program has 
been performed in the case where Ma ~ Ms ^ Mz and the results have been implemented 
in one of the RGE codes [20] . In the absence of such a tool for Ms ^ Ma ~ Mz (that is 
under development [21]), we will adopt the simple approach where the radiative corrections 
in the Higgs sector are approximated by the dominant contribution in the top and stop 
sector, which involves the logarithm of the scale Ms and the stop mixing parameter [22]. 
We will show that, in this approach, the situation simplifies to the extent that one can 
simply trade the dominant radiative correction against the actual value of the mass of the 
lighter h boson that has been measured at the LHC to be M^ ~ 126 GeV. 

One would then deal with a very simple post-/i discovery model in which, to a very 
good approximation, there are only two input parameters in the Higgs sector. Ma and 
tan/3 which can take any value (in particular low values tan/3 ~ 1 and Ma ~ 100 GeV 
unless they are excluded by the measurements of the h properties at the LHC) with the 
mass Mh fixed to its measured value. If one is mainly concerned with the MSSM Higgs 
sector, this allows to perform rather model-independent studies of this sector. 

We should note that while the working approximation for the radiative corrections to 
Higgs sector is important for the determination of the correct value of Ms (and eventually 
some other supersymmetric parameters such as the mixing in the stop sector), it has little 
impact on Higgs phenomenology, i.e. on the MSSM Higgs masses and couplings. 

The reopening of the low tan /3 region allows then to consider a plethora of very in- 
teresting Higgs channels to be investigated at the LHC: heavier CP-even H decays into 
massive gauge bosons H — )■ WW, ZZ and Higgs bosons H — )■ hh, CP-odd Higgs decays into 
a vector and a Higgs boson, A — )■ hZ, CP-even and CP-odd Higgs decays into top quarks, 
H/A — 7- ti, and even charged Higgs decay H — t- Wh. Many search channels discussed 
in the context of a heavy SM Higgs boson or for resonances in some non-SUSY beyond 
the SM (new gauge bosons or Kaluza-Klein excitations) can be used to search for these 
final states. A detailed discussion of the Higgs cross sections times decay rates in these 
process is made in this paper and an estimate of the sensitivity that could be achieved at 
the present -^i = 8 TeV run with the full data set is given. These processes allow to cover 
a large part of the parameter space of the MSSM Higgs sector in a model-independent 
way, i.e. without using the information on the scale Ms and more generally on the SUSY 
particle spectrum that appear in the radiative corrections. 

The rest of the paper is organised as follows. In the next section, we discuss the 
radiative corrections in the Higgs sector when Mh is used as input and their impact on the 
Higgs masses and couplings. In section 3, we summarize the various processes for Higgs 
production and decay in the high and low tan /3 regions and, in section 4, their implications 
for the MSSM parameter space. In section 5, we discuss the important new heavy Higgs 
channels that can be probed at the LHC at low tan/3. A conclusion is given in section 6. 



2. The Higgs sector of the MSSM in the various tan/3 regimes 

In this section, we review the theoretical aspects of the MSSM Higgs sector with some 
emphasis on the properties of the Higgs particles in the low tan/3 regime, 1 < tan/3 < 3, 
which contrary to the high tan /3 regime, has not received much attention in the literature. 

2.1 The radiatively corrected Higgs masses 

Let us begin by recalling a few basics facts about the MSSM and its extended Higgs sector. 
In the MSSM, two chiral superfields with respective hypercharges —1 and +1 are needed for 
the cancellation of chiral anomalies and their scalar components, the two doublet fields Hi 
and H2, give separately masses to the isospin — ^ and +^ fermions in a SUSY invariant way. 
After spontaneous symmetry breaking, the two doublet fields lead to five Higgs particles: 
two CP-even h,H bosons, a pseudoscalar A boson and two charged H^ bosons [3, 7]. 

The Higgs sector should be in principle described by the four Higgs boson masses and 
by two mixing angle a and /?, with a being the angle which diagonalises the mass matrix 
of the two CP-even neutral h and H states while /? is given in terms of the ratio of vacuum 
expectation values of the two Higgs fields Hi and H2, tan/3 = V2IV1. However, by virtue of 
SUSY, only two parameters are needed to describe the system at tree-level. It is common 
practice to chose the two basic inputs to be the pseudoscalar mass Ma-, expected to lie 
in the range between Mz and the SUSY breaking scale M5, and the ratio tan/?, which is 
expected to take values in the range [23] 

1 < tan/3<mi/7fi6«60 (2.1) 

with frit and m\, the running top and bottom quark masses in the MS renormalisation 
scheme evaluated at a scale close to the SUSY scale Ms- 

At tree-level, the CP-even h boson mass is then bound to be lighter than the Z boson, 
Mh < min(M^,Myi)| cos2/3| < Mz-, while the heavier H and H boson have masses that 
are comparable to that of the A state if M/i > Mz- Likewise, the mixing angle a can be 
written in compact form in terms of Mji and tan /3. If the mass Ma is large compared to 
the Z boson mass, the so called decoupling limit [11] that we will discuss in some detail 
here, the lighter h state reaches its maximal mass value, Mh ~ Mz\cos2l3\, the heavier 
CP-even and CP-odd and the charged Higgs states become almost degenerate in mass, 
Mh ~ Ma ~ Mh±, while the mixing angle a becomes close to a ~ | — /3. 

As is well known this simple pattern is spoiled when one includes the radiative cor- 
rections which have been shown to be extremely important [8,22,24-28]. Once these 
corrections are included, the Higgs masses (and their couplings) will, in principle, depend 
on all the MSSM parameters. In the phenomenological MSSM (pMSSM) [6], defined by 
the assumptions that all the soft-SUSY breaking parameters are real with the matrices 
that eventually describe them being diagonal (and thus, there is no new source of CP or 
flavor violation) and by the requirement of universal parameters for the first and second 
generation sfermions, the Higgs sector will depend on, besides Ma and tan/3, 20 additional 
parameters: the higgsino mass parameter n; the bino, wino and gluino mass parameters 
Mi,M2,M3; the first /second and third generation left- and right-handed sfermion mass 



parameters mq,muj^,m^ ^mj^men and mQ,mr,mj^ ,7?^j^,7?^f^; and finally the (common) 
first/second and third generation trilinear Au,Ad,Ae and At^Ab^A^- couplings^. 

Fortunately, only a small subset of these parameters has a significant impact on the 
radiative corrections to the Higgs sector. At the one loop level, the by far dominant 
correction to the Higgs masses is originating from top and stop loops and grows like the 
fourth power of the top quark mass, logarithmically with the stop masses and quadratically 
with the stop trilinear coupling. The leading component of this correction reads^ [22] 



3m^ 



4 



27r2t;2 sin^ /3 



1 M ^( _ ^t 



(2.2) 



where fht is again the running MS top quark mass to account for the leading two-loop QCD 
and electroweak corrections in a renormalisation group (RG) improvement (some higher 
order effects can also be included) [25]. We have defined the SUSY-breaking scale Ms to 
be the geometric average of the two stop masses Mg = ,/fnfTnf~; this scale is generally 
kept in the vicinity of the TeV scale to minimize the amount of fine tuning. We have also 
introduced the stop mixing parameter Xt = At — /icot/3, that we define here in the DR 
scheme, which plays an important role and maximizes the radiative correction when 

Xt = v6 Ms : maximal mixing scenario (2-3) 

while the radiative corrections is smallest for a vanishing Xt value, i.e. in the no mixing 
scenario Xt = 0. An intermediate scenario is when Xt is of the same order as the SUSY 
scale, Xt = Ms, the typical mixing scenario. These scenarios have been often used in the 
past as benchmarks for MSSM Higgs studies [31] and have been updated recently [32]. 

The e approximation above allows to write the masses of CP-even Higgs bosons in a 
particularly simple form 



MlH = l{Ml + Ml + e) 



/ M|M2 cos2 2/3 + e(M2 sin^ f3 + M^ cos2 (3) 



(2.4) 



(M2+M2+e) 

In this approximation, the charged Higgs mass does not receive radiative corrections, the 
leading contributions being of 0{amt) and one can still write the tree-level relation Mfj± = 
\/ M^ + M^. For large values of the pseudoscalar Higgs boson mass, the CP-even Higgs 
masses can be expanded in powers of M'^/M\ to obtain at first order 

2 „^„2 « ]\ ^-2 „■ 2 « I ^„^„2 
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(2.5) 



and indeed, in exact decoupling Ma/Mz —5- oo, one would have Mh = Ma= Mfj+ for the 
heavier Higgs states and, for the lighter h boson, the well known relation 



M/, = M;,"^^ = -^M|cos2 2/3+esin2/3 (2.6) 

^The first/second generation couplings have no impact in general and can be ignored in practice, reducing 
the effective number of free inputs of the pMSSM, from 22 to 19 parameters. 

^Note the typographical error for this equation in Ref. [7] which translated to Refs. [9, 29, 30]. 



In view of the large value M/j « 126 GeV of the observed Higgs state at the LHC, it is 
clear that some optimization of the various terms that enter the mass formula eq. (2.6) 
with the radiative correction eq. (2.2) is required. As was discussed in many instances 
including Refs. [9,29,30], one needs: i) to be close to the decoupling limit Ma ^ Mz 
and to have significant tan/3 values that lead to | cos 2/3 1 — t- 1 to maximize the tree-level 
mass and, ii) to be in the maximal mixing scenario Xt = y/GMg with the largest possible 
value of the SUSY-breaking scale Ms to maximize the radiative corrections. As the other 
SUSY-breaking parameters do not affect significantly the M™*^^ value, one can fix them 
to some value. For instance, one can make the choice [32] 

,. max . M2 :^ 2Mi = \fi\ = iMs, Ms = m,. = |m,-^ = I.5M5, A, = 0, 

where mq^ and m^. are the common first/second sfermion SUSY-breaking masses and Ai 
their trilinear couplings. Alternatively, one can perform a scan of these parameters in a 
reasonable range which should change the resulting value of MJ^^^ in the DR scheme only 
by a few GeV in general. 

In the case of a not too large SUSY scale, A/5 < 3 TeV, the numerical analyses of 
the MSSM Higgs sector can be performed with RGE programs [18, 19] such as Suspect 
which include the most relevant higher order radiative corrections in the calculation of the 
Higgs and superparticle masses (and their couplings). In particular, for the Higgs sector, 
the full set of one-loop radiative corrections which include also the sbottom and stau loop 
corrections that are important at high tan/3 values [24] and the dominant two-loop QCD 
and electroweak corrections [26] are incorporated in the DR scheme; the dominant three- 
loop corrections are also known [28] but they are quite small and they can be neglected. 

One should compare the results with those obtained with the program FeynHiggs [33] 
which incorporates the radiative corrections at the same level of accuracy but in the on- 
shell renormalisation scheme [27]. In most cases, one obtains comparable results but in 
some scenarios, the difference in the values of M/j can be as large as 3 GeV. We will thus 
assume, as in Ref. [32], that there is a AM^ ~ 3 GeV uncertainty on the determination of 
the h mass in the MSSM and that the value M^ = 126 GeV of the particle observed at the 
LHC corresponds to a calculated mass within the pMSSM of 

123 GeV <Mh< 129 GeV (2.8) 

This uncertainty includes the parametric uncertainties of the SM inputs, in particular the 
MS 6-quark mass and the top quark pole mass rnh{mi,) = 4:.7 GeV and m^° ° = 173.2ibl GeV 
[34]. In the latter case, it is assumed that the top quark mass measured at the Tevatron, 
with the uncertainty of 1 GeV, is indeed the pole mass. If the top mass is instead extracted 
from the top pair production cross section, which provides a theoretically less ambiguous 
determination of m^°^, the uncertainty would be of order 3 GeV [35]. Including also the 
experimental error in the Mh measurement by ATLAS and CMS, Mh = 125.7 it 0.4 GeV, 
the possible calculated mass value of the h boson in the MSSM can be extended to the 
much wider and admittedly rather conservative range 120 GeV<M/i<132 GeV. 



2.2 The low tan/3 regime 

The previous discussion assumed a not too high SUSY-breaking scale, Ms < 3 TeV, in 
order not to have a too large fine-tuning in the model. However, in many scenarios, 
values of Ms in the 10 TeV range and even beyond have been considered, with a most 
popular one being the split-SUSY scenario [14, 36]. Indeed, as the criterion to quantify the 
acceptable amount of tuning is rather subjective, one could well have a very large value of 
Ms which implies that no sfermion is accessible at the LHC or at any foreseen collider, with 
the immediate advantage of solving the flavor and CP problems in the MSSM by simply 
decoupling these states. The mass parameters for the spin-^ particles, the gauginos and 
the higgsinos, can be kept close to the electroweak scale, allowing for a solution to the dark 
matter problem and a successful gauge coupling unification, the two other SUSY virtues. 
The SUSY solutions to these two remaining problems are abandoned if one takes the very 
extreme attitude of assuming that the gauginos and higgsinos are also very heavy, with a 
mass close to the scale Ms-, as is the case of the so-called high-scale SUSY models [15, 36]. 

In all these these SUSY scenarios, there is still a light particle, the h boson, which can 
have a mass close to 126 GeV for a given choice of parameters such as Ms and tan/3; see 
for instance Refs. [9,36]. The other Higgs particles are much heavier as the pseudoscalar 
Higgs mass is very often related to the mass scale of the scalar fermions of the theory. 
Ma ~ Ms- However, this needs not to be the case in general, in particular for Ms values 
not orders of magnitude larger than 1 TeV. Even, in constrained minimal Supergravity- 
like scenarios, one can assume that the soft SUSY-breaking scalar mass terms are different 
for the sfermions and for the two Higgs doublets, the so-called non-universal Higgs mass 
models [16] in which the mass Ma is decoupled from Ms- Scenarios with very large values 
of Ms and values of Ma close to the weak scale have been advocated in the literature [17], 
while models in which one of the soft SUSY-breaking Higgs mass parameters, in general 
Mhi, is at the weak scale while Ms is large are popular; examples are the focus point 
scenario [37] and the possibility also occurs in M/string theory inspired scenarios [38]. 

Hence, if one is primarily concerned with the MSSM Higgs sector, one may be rather 
conservative and assume any value for the pseudoscalar Higgs mass Ma irrespective of 
the SUSY scale Ms- This is the quite "model-independent" approach that we advocate 
and will follow in this paper: we take Ma as a free parameter of the pMSSM, with values 
ranging from slightly above 100 GeV up to order Ms, but make no restriction on the SUSY 
scale which can be set to any value. 

Nevertheless, in scenarios with Ms ^ 1 TeV, the Higgs and SUSY mass spectrum 
cannot be calculated reliably using standard RGE programs as one has to properly decouple 
the heavy states from the low-energy theory and resum the large logarithmic corrections. 
A comprehensive study of the split SUSY spectrum has been performed in Ref. [20] and 
the various features implemented in an adapted version of the code SuSpect. However, 
this version does not include the possibility Ms ^ Ma ^ Mz that is of interest for us 
here. A comprehensive and accurate description of the high Ms scenario in the MSSM in 
the light of the h discovery, including the possibility of a Higgs sector at the weak scale, is 
under way [21]. In the meantime, we will use the e approximation of eq. (2.2) to describe 



the radiative corrections in our high Mg scenario which should be a good approximation 
for our purpose. In particular, for M4 ^ Mz, we have verified that our results are in a 
relatively good agreement with those derived in the more refined approach of Ref. [20] . 

Let us now discuss the magnitude of the SUSY scale that is needed to make small tan /3 
values viable. We make use of the program Suspect in which the possibility Ms ^ 1 TeV 
is implemented [20] and which includes the full set of radiative corrections (here we assume 
the maximal mixing Xt = ^/QMs scenario and we take 1 TeV for the gaugino and higgsino 
masses). In Fig. 1, displayed are the contours in the plane [tan/3, Ms] for fixed mass values 
Mh = 120, 123, 126, 129 and 132 GeV of the observed Higgs state (these include the 3 GeV 
theoretical uncertainty and also a 3 GeV uncertainty on the top quark mass). 
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Figure 1: Contours for fixed values Mh = 120, 123, 126, 129 and 132 GeV in the [tan/3, Ms] plane 
in the decoupling limit Ma ^ Mz] the "LEP2 contour" for Ady^ = 114 GeV is also shown. 

From the figure, one concludes that values of tan/3 close to unity are possible and allow 
for an acceptable M^ value provided the scale Ms is large enough. For instance, while one 
can accommodate a scale Ms ~ 1 TeV with tan/3 ~ 5, a large scale Ms ~ 20 TeV is 
required to reach tan/3 ~ 2; to reach the limit tan/3 = 1, an order of magnitude increase 
of Ms will be needed. Outside the decoupling regime, the obtained Ms for a given M^ 
value will be of course larger. For completeness, we also show the contour for the mass 
value Mh = 114 GeV, the 95% confidence level limit obtained at LEP2 on a SM-like Higgs 
boson; it illustrates the fact that values down to tan/3 ~ 1 are still allowed by this bound 
provided that Ms ^ 10 TeV. The implications of this feature will be discussed later. 

In the rest of this paper, we will thus consider situations with the MSSM Higgs sector 
at the weak scale and the only requirement that we impose is that it should be compatible 
with the LHC data and, in particular, with the mass and production rates of the Higgs 
boson that has been observed. The requirement that /Vf/i~ 126 GeV, within the theoretical 
and experimental uncertainties, will be turned into a requirement on the parameters that 
enter the radiative corrections and, hence, on the scale Ms and the mixing parameter Xt, 
for given values of the two basics inputs Ma and tan/3. 



2.3 The Higgs couplings and the approach to the decouphng hmit 

Let us now turn to the important issue of the Higgs couplings to fermions and gauge 
bosons. These couplings strongly depend on tan/3 as well as on the angle a (and hence on 
Ma)', normalized to the SM Higgs couplings, they are given in Table 1. The A boson has 
no tree level couplings to gauge bosons as a result of CP-invariance, and its couplings to 
down-type and up-type fermions are, respectively, proportional and inversely proportional 
to tan/3. This is also the case for the couplings of the charged Higgs boson to fermions, 
which are admixtures of 7fifetan/3 and ifit cot j3 terms and depend only on tan/3. For the 
CP-even Higgs bosons h and H, the couplings to fermions are ratios of sines and cosines of 
the angles a and /?; the couplings to down (up) type are enhanced (suppressed) compared 
to the SM Higgs couplings for tan/3 > 1. The two states share the SM Higgs couplings 
to vector bosons as they are suppressed by sin(/3 — a) and cos(/3 — a), respectively for 
h and H. We note that there are also couplings between a gauge and two Higgs bosons 
which in the case of the CP-even states are complementary to those to two gauge bosons 
QhAZ oc QhH+w- ^ dHVV and vice versa for h -^ H; the coupling gAH+w- has full strength. 
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Table 1: The couplings of the neutral MSSM Higgs bosons, collectively denoted by <&, to fermions 
and gauge bosons when normalized to the SM Higgs boson couplings. 

These couplings are renormalized essentially by the same radiative corrections that 
affect the CP-even neutral Higgs masses. In the e approximation discussed above, the 
one-loop radiatively corrected mixing angle a will indeed read 



tan 2a = tan 2/3 



Ml + Ml 



IT 
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Ml + e/ cos 2/3 
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(2.9) 



This leads to corrected reduced h, H couplings to gauge bosons that are simply ghvv = 
sin(/3 — a) and gnvv = cos(/3 — a) and similarly for the couplings to fermions. 

The decoupling limit is controlled by the VV coupling of the heavier CP-even Higgs 
boson, gHVV = cos(/3 — a), which vanishes in this case, while the hVV coupling gf^yy = 
1 — g^yy = sin (/3 — a) becomes SM-like. Performing again an expansion in terms of the 
pseudoscalar Higgs mass, one obtains in the approach to the decoupling limit^ 
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(2.10) 



where, in the intermediate step, the first term is due to the tree-level contribution and the 
second one to the one-loop contribution e. Concentrating first on the tree-level part, one 



^We thank Nazila Mahmoudi for discussions concerning these Umits. 



realises that for large values of tan /? and also for values very close to unity, the decoupling 
limit is reached more quickly. Indeed the expansion parameter involves also the factor 
sin 4/3 which becomes in these two limiting cases 

. ,^ 4tan/3(l -tan^m f -4/ tan/3 for tan/3 >1 ^ ,^^^, 

sm4/3 = — -— ^ „ ^,^ -» < o -^0 (2.11) 

^ (l + tan2/3)2 \l-tan2/3for tan/3 ~ 1 ^ ^ 

Hence, in both the tan/3 ^ 1 and tan/3 ~ 1 cases, the gnw coupling that controls the 
decoupling limit M'^/M\ — )• 0, is doubly suppressed. The radiatively generated component, 
if one recalls that the one-loop correction in eq. (2.2) involves a 1/ sin^ j3 term which makes 
it behave as —e/M\ x cot /3, also vanishes at high tan /3 values. This leads to the well known 
fact that the decoupling limit gnw — )• is reached very quickly in this case, in fact as soon 
as Ma ^ M™^^. Instead, for tan/3 ~ 1, this radiatively generated component is maximal. 
However, when both components are included, the departure from the decoupling limit in 
the coupling gnw for a fixed Ma_ value occurs when sin 4/3 ~ —1, which corresponds to 
/3 = 37r/8 and hence to the value tan/3 ss 2.4. 

Similarly to the HVV case, one can write the couplings of the CP~even Higgs states 
to isospin ^ and —^ fermions in the approach to the decoupling limit M^/M^ <C 1 as 

ghuu — > 1 + X cot /3 -^l 

Ma->Mz . J. « . 1 

ghdd — > 1 - X tan /3 ^ 1 

gHuu — > - cot /3 + X -^ - cot /3 

gHdd — > +tan^ + x -;>+tan^ 

with the expansion parameter x ^ ^/^\ is the same as the one given in eq. (2.10). In the 
Ma ^ Mz regime, the couplings of the h boson approach those of the SM Higgs boson, 
9huu ~ ghdd ~ 1) while the couplings of the H boson reduce, up to a sign, to those of 
the pseudoscalar Higgs boson, guuu » gAuu = cot/3 and gudd ~ gAdd = tan/3. Again, as 
a result of the presence of the same combination of M^ sin 4/3 and e sin 2/3 factors in the 
expansion term x of all couplings, the limiting values are reached more quickly at large 
values of tan/3 but the departure from these values is slower at low tan/3. 

In Fig. 2, we display the square of the H couplings to gauge bosons and fermions as a 
function of tan/3 for Ma = 300 GeV. Again the maximal mixing scenario is assumed and 
Ms is chosen in such way that for any tan/3 value, one has M^ = 126 GeV. At such A 
masses, the couplings of the lighter h boson to all particles deviate little from unity even 
for small tan/3 values and in this case too one can consider that we are already in the 
decoupling regime. Nevertheless, the coupling of the heavier H boson to VV states is still 
non-zero, in particular at low tan/3. The H coupling to ti pairs states (as well as the A 
coupling) is significant at low tan/3 values, g'^^^ > 0.1 for tan/3 < 3. 

This demonstrates that the heavier H/A/H bosons can have sizable couplings to top 
quarks (and to massive gauge bosons for H outside the decoupling regime) if tan /3 values 
as low as tan/3 ~ 3 are allowed. In fact, the H/A/H^ couplings to top quarks oc cot/3 are 
larger that the couplings to bottom quarks oc tan /3 for values tan /3 w yjfht/rhi, < 7 and this 
value should be considered as the boundary between the high and low tan/3 regimes. With 
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Figure 2: The squared couplings of the heavier CP-even H state to gauge bosons and fermions as 
a function of tan;3 for Ma = 300 GeV. The SUSY scale is chosen so that Mn = 126 GeV. 

more refinement, one can consider three tan/3 regimes: the high regime with tan/3 > 10, 
the intermediate regime with 5 < tan/3 < 10 and the low regime with tan/3 < 5. 

There are two important remarks which should be made before closing this section. The 
first one is that besides the e correction, there are additional one-loop vertex corrections 
which modify the tree-level Higgs-fermion couplings [39]. In the case of 5-quarks in the 
high (and eventually intermediate) tan/3 regime, they can be very large in the 6-quark 
case as they grow as mfttan/3. The dominant component comes from the SUSY-QCD 
corrections with sbottom-gluino loops that can be approximated by 

Afo ~ — -firrig tan /3/Taax{m^,mr , m,r ) (2-12) 

In the decoupling limit Ma ^ Mz , the reduced bb couplings of the H, A states read 

1 + A6_ 

In the case of the lighter h boson, the hbb couplings stay SM-like in this limit in principle, 
but slightly outside the decoupling limit, there is a combination of SUSY parameters which 
realises the so-called "vanishing coupling" regime [31] in which a — )• and hence g^bb "^ 1- 
The second remark concerns the trilinear Hhh coupling which will be needed in our 
analysis. In units of M'^/v, this coupling is given at tree-level by [3] 

^Hhh ~ 2 sin 2a sin(/3 + a) — cos 2q cos(/3 + a) (2-14) 

Again, to include the radiative corrections in the e approximation, one needs to perform 
the change a — t- a; however, in this case, there are also direct vertex corrections but they 
can be still described by the e parameter. One obtains in this approach [40] 

3 



gHbb ~ gAbb ~ tan /3 



1 



(2.13) 



^Hhh 



Ma:^Mz 



Ml 



{Ml-e sin^ /3) ( M| - M^ + e sin^ /3) + e sin /3 cos /3 



(2.15) 



At high-tan /3, the trilinear coupling vanishes Xuhh ~^ while for small and intermediate 
tan /3 values it stays quite substantial as a result of the large e corrections. 
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3. Higgs decays and production at the LHC 

3.1 The high and intermediate tan/3 regimes 

The production and decay pattern of the MSSM Higgs bosons crucially depend on tan/3. 
In the LHC run up to now, i.e. with center of mass energies up to \/s = 8 TeV, only 
relatively large tan/3 values, tan/3 > 5-10 which correspond to the high and intermediate 
regimes, are probed in the search of the neutral H/A and the charged H bosons. In the 
high tan/3 regime, the couplings of these non-SM like Higgs bosons to h quarks and to 
r leptons are so strongly enhanced, and the couplings to top quarks and massive gauge 
bosons so suppressed, that the pattern becomes rather simple. 

A first simplifying feature is that the decoupling regime in which the lighter h boson 
attains its maximal mass M™^^ value for a given SUSY parameter set^ and has SM- 
couplings is reached already at Ma ^ M™**^ for tan/3 > 10. In this case, the heavier CP- 
even H boson has approximately the same mass as the A boson and its interactions are 
similar. Hence, the spectrum will consist of a SM-like Higgs h = H^m and two pseudoscalar 
(like) Higgs particles, <I> = H/A. The H boson will also be approximately degenerate in 
mass with the <I> states and the intensity of its couplings to fermions will be similar. 

An immediate consequence will be that the h boson will precisely decay into the variety 
of final states and will be produced in the various channels that are present in the SM. 
These decay and production processes have been studied in detail at various places, see 
Ref. [4] for a detailed review and Refs. [42, 44] for updates. We will discuss the implications 
of these channels for the properties of the state observed at the LHC in the next section. 

In the case of the heavier neutral <I> = H/A bosons, the decay pattern is very simple: 
the ti channel and all other decay modes are suppressed to a level where their branching 
ratios are negligible and the <I> states decay almost exclusively into t~^t~ and hb pairs, with 
branching ratios of BR($ -^ t'^t') « ml/[?,ml{Mq>)+ml] « 10% and BR($ -^ hi) ^ 90%. 
The charged Higgs particles decay into H — t- rvr final states with a branching fraction of 
almost 100% for H^ masses below the tb threshold, Mjj± > mt + nih, and a branching ratio 
of only ~ 10% for masses above this threshold. The by far dominant channel in the latter 
case is H^ — )■ tb which occurs with a ~ 90% probability for the same reason as above. 

Concerning Higgs production in the high tan/3 regime, the enhancement of the 6-quark 
couplings makes that only processes involving this quark are important for the <I> = H/A 
states. In the dominant gluon fusion production channel, gg — t- <5, one should take into 
account the 6-quark loop which provides the largest contribution (in contrast to the SM 
where the top quark contribution largely dominates) and in associated Higgs production 
with heavy quarks, bb final states and hence the processes gg/qq — )• bb + ^, must be 

*The present discussion holds in the case where the h boson is the SM-like state which implies 
Ma > M™'^^. At low Ma values, the role of the CP-even h and H states are reversed: it is H which 
is the SM-like particle H = Hsm and h would correspond to the pseudoscalar-like Higgs particle. However, 
the possibility that the H state is the observed particle at the LHC is ruled out by present data [30]. A 
special case would be Ma ~ M™^", which is called the intense coupling regime in Ref. [41] and which leads 
to mass degenerate h, H, A states with comparable couplings to fermions; as the h and H states are close 
in mass, one has the same phenomenology as in the decoupling limit where H has the same properties as 
A [42,43]. Again, this scenario is strongly disfavored by present data [30]. 



12 



considered. The latter processes are equivalent to the 66 —t- <^ channels when no-additional 
6-quark in the final state is present, if one considers the 6-quark as a massless parton and 
uses heavy quark distribution functions in a five active flavor scheme [45]. 

Hence, except for the gg ^ ^ and 66 — )■ <& fusion processes, all the other production 
channels are irrelevant in the high tan/3 regime, in particular the vector boson fusion and 
the Higgs-strahlung channels, that are absent for A and strongly suppressed for H. In 
both cases, as M$ ^ rrih^ chiral symmetry holds and the cross sections are approximately 
the same for the CP-even H and CP-odd A bosons. The cross section for gg ^^ ^ \s 
known up to next-to-leading order in QCD [46] and can be calculated using the program 
HIGLU [47,48]. The 66 ^ $ rate is instead known up to NNLO in QCD [49,50] and its 
evaluation can be made using the programs bbSnnlo or SUSHI [51]. Note that for associated 
H/A production with two tagged 6-quarks in the final states that can be used, one should 
instead consider the process gg/qq — )■ 66 + $ which is known up to NLO QCD [52]; they 
leading order cross section can be obtained using the program QQH [48]. 

The most powerful search channel for the heavier MSSM Higgs particles at the LHC 
is by far the process 

pp -^ gg + bb -?■ ^ -?- T^T^ (3.1) 

The precise values of the cross section times branching fraction for this process at the LHC 
have been recently updated in Refs. [42, 44] and an assessment of the associated theoretical 
uncertainties has been made. It turns out that these uncertainties are not that small. They 
consist mainly of the scale uncertainties due to the missing higher orders in perturbation 
theory and of the combined uncertainty from the parton distribution functions and the 
strong coupling constant as- When combined, they lead to a total theoretical uncertainty 
of 20-30% in both the gg ^ ^ and 66 — )■ <I> channels^. We will assume here for the combined 
gg + bb ^ ^ channel a theoretical uncertainty of 

A™o-(pp^$) X BR($^ rr) = ±25% (3.2) 

in the entire M$ range probed at the LHC and for both i/s = 8 and 14 TeV. 

Besides the QCD uncertainty, three other features could alter the rate a{pp^-^^TT) 
in the MSSM and they are related to the impact of the SUSY particle contributions. We 
briefly summarise them below and some discussions are also given in Refs. [30,53]. 

While the CP-odd A state does not couple to identical squarks as a result of CP- 
invariance, there is a Hqiqi coupling in the case of the H state which allows squarks, 
and mainly top and bottom squarks, to contribute to the gg ^ H amplitude at leading 
order (there are NLO contributions [54] for both the Hgg and Agg amplitudes via gluino 
exchange but they should be smaller). However, as squarks do not couple to the Higgs 
bosons proportionally to their masses, these contributions are damped by powers of m^ 



^It was advocated in Refs. [42, 43] that there are two additional sources of uncertainties related to 
the 6-quark mass which should be considered: the one in the gg -^ ^ process due to the choice of the 
renormalization scheme for nib and the parametric uncertainty. These will increase the total uncertainty. 
We will however, ignore this complication and retain the "official" estimate of the error given in Ref. [44]. 
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for Mu < 2mQ and, at high tan/3, the 6-loop contribution stays largely dominant. These 
SUSY contributions are thus expected to be small and can be neglected in most cases. 

A more important effect of the SUSY sector is due to the one-loop vertex correction 
to the ^bb couplings, Af, of eqs. (2.12-2.13), which can be large in the high tan/3 regime 
as discussed previously. However, in the case of the full process pp — ?• $ — )■ t~^t~ , this 
correction appears in both the cross section, (t(<1>) oc (1 + Af,)"^, and in the branching 
fraction, BR(rr) = r($ -)■ tt)/[(1 + Afe)-2r($ -^ bb) + r($ -^ tt)], which involves the 
Ab correction above in the denominator. Hence, in the cross section times branching ratio, 
the Ab corrections largely cancels out and for BR(tt) ks 10%, one obtains 

(T{gg + bb^<^)x BR($ -^ tt) ^ a x BR x {1 - -Afe) (3.3) 

o 

Hence, one needs a very large A;, term (which, one should recall, is a radiative correc- 
tion), of order unity or more, in order to alter significantly the pp — )•<!>—)• rr rate^. 

Finally, there is the possibility that there are light SUSY particles with masses m < ^M^ 
which lead to the opening of SUSY decay channels for the H/A states that might reduce 
the $ — 7- rr branching fraction. For M$ < 1 TeV, the only possibilities for these superpar- 
ticles seem to be light neutralinos and charginos (x) and light sleptons (£). These decays 
have been reviewed in Ref. [7] and they have been found to be in general disfavored in 
the high tan (3 regime as the $ — t- 66 + rr decays are so strongly enhanced that they leave 
little room for other possibilities. Only in a few special situations that these SUSY decays 
can be significant. For the decays ^ — t- xX^ it is the case when i) all x = Xi 2 ^^"^ /ti-4 
channels are kinematically open or ii) if only a subset of x particles is light, they should 
be mixtures of gauginos and higgsinos to maximize the $XX couplings. Both scenarios 
should be challenged by the present LHC constraints^. In the case of sleptons, only the 
decays into light f states could be important; while the decay A — )• fifi is forbidden by 
CP-invariance, the decays H — )■ fifi and H/A — )• fif2 can have substantial rates at high 
tan/3 when the '^ff coupling is enhanced. However, again, at these large tan/3 values, the 
$ —7- 66 and ^ ^ tt decays are extremely enhanced and leave little room for competition. 

Thus, only in the unlikely cases where the decay ff — >-'fifi has a branching rate of the 
order of 50%, the squark loop contribution to the gg ^ H process is of the order 50%, or 
the Afo SUSY correction is larger than 100%, that one can change the pp — )■ $ — )■ rr rate by 
~ 25%, which is the level of the QCD uncertainty. One thus expects cr(pp — )• <I>) x BR($ — )• 
tt) to be extremely robust and to depend almost exclusively on Ma and tan/3. 

Two more processes are considered for the heavier MSSM neutral Higgs bosons at 
high tan/3. The first one is pp — t- <I> — t- /i^/i~ for which the rate is simply cr(pp— t-'I' — t- tt) 



"in any case, if one insists to take this A;, correction into account in tlie constraint on the [tan /3, Ma] 
plane that is obtained from the pp —>■$—>■ rr rate, one could simply replace tan/? by tan/J/(l + Ai,/10). 
A contribution A^ ~ 1 will change the limit on tan/? by only 10%, i.e. less than the QCD uncertainty. 

^The searches of charginos and neutralinos in the same-sign lepton and tri-lepton topologies at the LHC 
are now probing significant portions of the gaugino-higgsino parameter space and they exclude more and 
more the possibility of light x states [55]. This is particularly true for mixed gaugino-higgsino states in 
which the <i>xx couplings are maximised: the lead to a large gap between the lightest and the next-to-lightest 
X masses and hence a large amount of missing energy that make the searches more effective. 
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rescaled by BR(<I> — t- ;U/x)/BR(<I> — t- tt) = m? /m^ ~ 4 x 10~^. The rate is much sinaUer 
than in the tt case and is not compensated by the much cleaner ^^ final state and the 
better resolution on the invariant mass. Searches in this channel have been performed in 
the SM Higgs case [56] and the sensitivity is very low. In addition, there is the process in 
which the H/A bosons are produced in association with two 6-quark jets and decay into 
hh final states and searches in this channel have been performed by the CMS collaboration 
with the 7 TeV data [57]. However, the sensitivity is far lower than in the t~^t~ channel. 

Thus, the pp ^ $ — )■ t~^t~ search for the neutral Higgs bosons provides the most 
stringent limits on the MSSM parameter space at large tan/3 and all other channels are 
weaker in comparison and provide only cross checks. We will thus concentrate on this 
process in the rest of our discussion of the high tan /3 regime. 

A final remark needs to be made on the charged Higgs boson. The dominant H^ search 
channel at present energies is in H^ — )• tu final states with the H^ bosons produced in 
top quark decays for masses not too close to Mfj± = mt — mf, ss 170 GeV 

pp -^ tt with t -^ H+h ^Tvh (3.4) 

This is particularly true at high tan/3 values when the t — )■ H^b branching ratio which 
grows with m^ tan^ /3, is significant. For higher H masses, one should rely on the three- 
body production process pp — )• tbH — )• tbTv which leads to a cross section that is also 
proportional to tan^ /3, but the rates are presently too small. Hence, processes beyond 
t — )• bH^ can be considered only at the upgraded LHC. 

3.2 The lovif tan/3 regime 

The phenomenology of the heavy MSSM A, H, H^ bosons is richer at low tan /3 and leads 
to a production and decay pattern that is slightly more involved than in the high tan/3 
regime. Starting with the production cross sections, we display in Fig. 3 the rates for the 
relevant H/A/H production processes at the LHC with center of mass energies of y/s = 8 
TeV and ^ = 14 GeV assuming tan/3 = 2.5. The programs HIGLU [47], SUSHI [51] and 
those of Ref. [48] have been used with a calling file in which the radiative corrections in 
the Higgs sector are calculated according to section 2.1 and lead to a fixed M^ = 126 GeV 
value. The MSTW set of parton distribution functions (PDFs) [58] has been adopted. For 
smaller tan/3 values, the cross sections for the various processes, except for pp — t- H/A + bb, 
are even larger as the H/A couplings to top quarks and the HVV coupling outside the 
decoupling limit are less suppressed. 

Because of CP invariance which forbids AVV couplings, the pseudoscalar state A can- 
not be produced in the Higgs-strahlung and vector boson fusion processes. For Ma ^ 300 
GeV, the rate for the associated pp — )• tiA process is rather small, as is also the case of 
the pp —7- bbA cross section which is not sufficiently enhanced by the Abb oc tan /3 coupling. 
Hence, only the gg ^ A fusion process with the dominant t-quark and sub-dominant 
6-quark loop contributions included provides large rates at low tan/3. 

The situation is approximately the same for the CP-even H boson: only the gg ^ H 
process provides significant production rates at relatively high values of Mh, Mh > 300 
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GeV, and low tan/3, tan/3 < 5. As in the case of A, the cross section for pp — t- ttH is 
suppressed compared to the SM case while the rate for pp — )• bbH is not enough enhanced. 
However, in this case, the vector boson fusion pp — t- Hqq and Higgs-strahlung processes 
qq —7- HW/HZ are also at work and have production rates that are not too suppressed 
compared to the SM at sufficiently low A4h values, Mh < 200-300 GeV and tan/3 ~ 1. 
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Figure 3: The production cross sections of the MSSM heavier Higgs bosons at the LHC with 
^/s — 8 TeV (left) and ^ = 14 TeV (right) for tan/3 = 2.5. Only the main channels are presented. 
The higher order corrections are included (see text) and the MSTW PDFs have been adopted. 

Hence, for M^ ^ 300 GeV, the only relevant production process is gg ^ ^ with the 
dominant contribution provided by the heavy top quark loop. In this case, one can include 
not only the large NLO QCD corrections [59], which are known in the exact case [46], but 
also the NNLO QCD corrections [60] calculated in a effective approach with mt ^ M$ 
which should work in practice for M$ < 300 GeV but can be extended to higher masses. 

For the charged Higgs boson, the dominant production channel in the low tan/3 regime 
is again top quark decays, t — )■ H~^b, for Mfj± < 170 GeV. Indeed, for tan/3 < 5, the 
mt/ tan /3 component of the H tb coupling becomes rather large, leading to a significant 
t —7- H^b branching ratio. For higher H^ masses, the main process to be considered is 
99 /qq ~^ H^tb [61]. As in the case oi pp — )• bb^, one can take the 6-quark as a parton and 
consider the equivalent but simpler 2 — )• 2 channel gb — >• H^t. One obtains an accurate 
description of the cross section if the renormalisation and factorisation scales are chosen to 
be low, fiji = fip ^ ^{Mij± + mt) in order to account for the large NLO QCD corrections 
[62]; the scales uncertainties are large though, being of order 20% [44]. Additional sources 
of H^ states for Mh± < 250 GeV are provided by pair and associated production with 
neutral Higgs bosons in qq annihilation as well as H'^H" pair and associated H W^ 
production in gg and/or bb fusion but the rates are very small [63] . 

Let us turn to the decay pattern of the heavier MSSM Higgs particles which can be 
rather involved in the low tan/3 regime. In this case, as the couplings of the H/A bosons 
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Figure 4: The decay branching ratios of the heavier MSSM Higgs bosons A (left), H (center) and 
H^ (right) as a function of their masses for tan/3 — 2.5. The program HDECAY [64] has been used 
with modifications so that the radiative corrections lead to Mh = 126 GeV. 

to 5-quarks are not very strongly enhanced and the couplings to top quarks (and gauge 
bosons in the case of the H state) not too suppressed, many interesting channels appear. 
The branching fractions for the H/A/H^ decays are shown in Fig. 4 as functions of their 
masses at tan/3 = 2.5. They have been obtained using the program HDECAY [64] assuming 
large Mg values and the maximal mixing scenario that lead to a fixed M^ = 126 GeV value. 
The pattern does not significantly depend on other SUSY parameters, provided that Higgs 
decays into supersymmetric particles are kinematically closed as will be implicitly assumed 
in the following^, where the main features of the decays are summarised in a few points. 

- Sufficiently above the ti threshold for the neutral and the th threshold for the charged 
Higgs bosons, the decay channels H/A — )• ti and H^ — )• th become by far dominant for 
tan/3 < 3 and do not leave space for any other decay mode. Note that these decays have 
also significant branching fractions below the respective kinematical thresholds [65]. It is 
especially true for the charged Higgs state for which BR(ff + — )• t6) > 1% for Mh± w 130 GeV. 

- Below the ti threshold, the H boson can still decay into gauge bosons H — t- WW and 
ZZ with rather substantial rates as the HVV couplings are not completely suppressed. 

- In the window 2M/j < Mh < 2mt, the dominant decay mode for tan/3 < 3 turns 
out to be the very interesting channel H ^ hh channel. As discussed earlier, the Hhh 
self-couplings given in eq. (2.15) is significant at low tan/3 values. 

- If allowed kinematically, i.e. for AIa ^ ^Ih + Mz GeV, the CP-odd Higgs boson can 
also decay into hZ final states with a significant rate below the ti threshold as the AZh 
coupling (that is the same as the HVV coupling) is substantial. Nevertheless, the A ^ tt 
channel is still important as it has a branching fraction above ~ 5% up to Ma ~ 2mt. 



*In fact, even in this low tan /J case, the ti decays for sufficiently large masses are so dominant that 
they do not lead to any significant quantitative change if SUSY particles are light. In addition, being not 
enhanced by tan/3, the A;, correction has no impact in this low tan/3 regime. 
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- In the case of the charged Higgs state, there is also the channel H^ — )• Wh which is 
important similarly to the A — >• hZ case. Note that for Mu± < 170 GeV, the decay H^ — )■ 
cs that is usually considered only in two-Higgs doublet models and the very interesting 
flavor changing mode H'^ — )• cb have rates that are at the percent level. All these exotic 
channels have larger branching ratios, above ~ 10%, for tan/3 values close to unity. 

3.3 The case of the h boson 

A final word is due to the production and decay rates of the lighter h boson that we will 
assume to be the state with a 126 GeV mass observed at the LHC. 

In the left-hand side of Fig. 5, we display the cross sections for the relevant Higgs 
production channels at the LHC with -^s = 8 TeV as a function of M^ at tan/3 = 2.5. 
Again, the radiative corrections in the e approach are such that M^ is fixed to 126 GeV. 
Shown are the rates for the gluon fusion gg ^ h, vector boson fusion qq — )• Hqq, Higgs- 
strahlung qq — )• hW, hZ as well as associated pp — t- tth processes. The relevant higher order 
QCD corrections are implemented and the MSTW set of PDFs has been adopted. The 
rates can be very different whether one is in the decoupling limit M^ ~ 1 TeV where the 
h couplings are SM-like or at low M^ values when the h couplings are modified. 

The variation of the branching ratios compared to their SM values, which correspond to 
their IVISSM values in the decoupling limit, are displayed as a function of Ma_ for tan /3 = 2.5 
in the right-hand side of the figure. Sown are the branching fractions for the decays that are 
currently used to search for the SIVI Higgs boson, i.e. the channels h — )• bb, tt, ZZ, WW, 77. 
Again, large differences compared to the SM can occur at low to moderate Mj^ values. 









1.8 

1.6 

i 1.4 
B 
S" 1.2 

ffl 
J 1 

CQ 0.8 
0.6 
0.4 


h^bb 

tan^ = 2.5 h^rf 

h^77 

. Mh = 126GeV h^ZZ 




10 








ft 1 


h^WW 


J3 


"~^^ 


t 

ft °'^ 

0.01 
0.001 


' 










— ^ — 




V^ = 8TcV 
- tan/3 = 2.5 
Mh = 126 GeV 


ggh 

hqq 

bbh 
Wh 
Zh 
tth 


- 


^^'■"^ 


^ 


200 


400 600 
Ma [GeV] 


m ic 


00 


300 400 600 800 
Ma [GeV] 


10 



Figure 5: The production cross sections of the lighter h boson at the LHC with -y/s = 8 TeV (left) 
and the variation of its decay branching fractions compared to the SM values (right) for tan/3 = 2.5. 
Again, the radiative corrections in the Higgs sector are such that Af/j = 126 GeV. 

The data collected so far by the ATLAS and CMS collaborations on the observed 126 
GeV Higgs particle should thus put strong constraints on the parameters tan/3 and Afy^. 
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4. Present constraints on the MSSM parameter space 

4.1 Constraints from the h boson mass and rates 

We start this section by discussing the impact of the large amount of ATLAS and CMS 
data for the observed Higgs state at the LHC on the MSSM parameter space. We will 
assume for definiteness that the h boson is indeed the observed particle as the possibility 
that it is the H state instead is ruled out by the LHC data [30]. 

A first constraint comes from the measured mass of the observed state, Mh ~ 126 
GeV. As discussed previously and in several other instances such as Ref. [9], in the phe- 
nomenological MSSM, this large Mh value indicates that the radiative corrections in the 
Higgs sector are maximised. If the scale Ms is close to 1 TeV as dictated by naturalness 
arguments, this implies that one is in the decoupling regime (and hence, dealing with a 
SM-like Higgs particle) with intermediate to high-tan /3 values and maximal stop mixing. 
If the SUSY scale is pushed to Ms ~ 3 TeV, the highest acceptable value from fine-tuning 
adopted in many analyses such that of Refs. [29, 30], a smaller mixing in the Higgs sector 
and values of M^ of order of a few hundred GeV can be made possible. However, tan/3 
values in the low regime, tan/3 < 3-5 cannot be accommodated as they lead to Mh < 123 
GeV and even to Mh < 120 GeV, which is the lowest value that can be reached when 
including the theoretical and the top-quark mass uncertainties in the calculation of Mh- 

To obtain an acceptable value of Mh in the low tan /3 regime, one needs to push Ms to 
the 10 TeV domain or higher. In the approach that we are advocating here, in which the 
radiative corrections in the MSSM Higgs sector are implemented in the rather simple (but 
not completely inaccurate) approximation where only the leading RGE improved one-loop 
correction of eq. (2.2) is taken into account, one can trade the (unknown) values of Ms and 
the mixing parameter Xt with the (known) value of the Higgs mass Mh- In other words, 
for each set of tan/3 and M^ inputs, one selects the e radiative correction that leads to the 
correct mass Mh = 126 GeV. The LHC constraint on the mass of the observed Higgs state 
is then automatically satisfied. We emphasize again that for the large SUSY scales that 
are needed for the low tan/3 regime, tan/3 < 3, the MSSM spectrum cannot be calculated 
in a reliable way using the usual versions of the RGE programs such as Suspect. 

A second constraint comes from the measurement of the production and decay rates 
of the observed Higgs particle. The ATLAS and CMS collaborations have provided the 
signal strength modifiers ^xx, that are identified with the Higgs cross section times decay 
branching ratio normalized to the SM expectation in a given H^XX search channel. For 
the various search channels that have been conducted, h — )■ ZZ, WW, 77, tt and bb with 
the entire set of data collected in the runs at ^/s = 7 TeV and 8 TeV, i.e. ~ 5 fb~^ and 
~ 20 fb~^ (with the exception oi h ^ bb which has been analyzed only with 17 fb~^ of the 
7-1-8 TeV data) numbers can be found in Refs. [66-75]. These should be used to constrain 
the couplings of the h state and, hence, the [tan/3,Myi] parameter space. 

Rather than performing a complete fit of the ATLAS and CMS light Higgs data includ- 
ing all the signal strengths, we will simply use the most precise and cleanest observable in 
this context: the signal strength nzz in the search channel h — ?■ ZZ. As recently discussed 
in Ref. [78] (to which we refer for the details), this channel is fully inclusive and does not 
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involve the additional large theoretical uncertainties that occur when breaking the cross 
section of the dominant production process gg^h into jet categories^. In addition, con- 
trary to the global signal strength //tot, it does not involve the channel /i — t- 77 which, at 
least in the ATLAS case, deviates from the SM prediction and might indicate the presence 
of new contributions (such as those of light charginos?) in the h'j'j loop. The combination 
of the ATLAS and CMS data in the ZZ channel gives, //zz = 1, 10 ± 0.22 ± 0.20 where the 
first uncertainty is experimental and the second one theoretical. Following Ref. [42], we 
assume a total theoretical uncertainty of A*^ = ±20% and, since it should be considered as 
a bias, we add it linearly to the experimental error. This gives a lower limit on the h — )> ZZ 
signal strength of fizz > 0.62 at 68%CL and fizz > 0.4 at 95%CL. 

In the MSSM case, the signal strength will be given by fizz = c(/i) x BR(/i — )• ZZ)/ 
o'iHsm) X BR(i?sM — ^ ZZ) and will be thus proportional to combinations of reduced h 
coupling squared to fermions and gauge bosons, g^^j x ghvv/dhbb--- "^^^ ^^^^ ^^^* A*^-^ ^^^^ 
be as low at 0.4 at 95%CL means that we can be substantially far from the decoupling 
limit, gfjyy ~ 0.1, with not too heavy H/A/H^ states even at low tan/3. 

In Fig. 6, we have scanned the [tan/3,M^] parameter space and delineated the areas 
in which the 68%CL and 95%CL constraints on fizz are fulfilled. We observe that indeed, 
the entire range with Ma ^ 200 GeV for most value of tan/3 is excluded at the 95%CL. 
With increasing tan/3, the excluded M4 values are lower and one recovers the well known 
fact that the decoupling limit is reached more quickly at higher tan /3 values. In most cases, 
we will use this indirect limit of Ma ^ 200 GeV prior to any other constraint (except for 
illustrations in the H^ case where the low mass range will be kept). 
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Figure 6: The [tan/3, M^] parameter space of the MSSM in which the signal strength in the 
h -^ ZZ search channel is compatible with the LHC data on the rates of the observed h boson at 
the 68%CL (green), 95%CL (yellow) and 99%CL (blue). 



®For instance, the signal strengths is the tt and WW channels are obtained by considering the gg — )> H-\- 
Oj, Ij and/or the vector boson fusion categories. The signal strength ^ww provides the same information 
as fizz, while the measurement of the signal strengths in the h ^ bb and h — >■ t^t^ channels are not yet 
very accurate. Hence, using only the h — >■ ZZ channel should be a good approximation. 
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4.2 Constraints from the heavier Higgs searches at high tan/3 

As discussed in section 3.1, tlie most efficient channel to search for the heavier MSSM 
Higgs bosons is by far H/A production in gg and hh fusion with the Higgs bosons decaying 
into r lepton pairs, pp — t- <I> — )■ t^t^ . Searches for this process have been performed 
by the ATLAS cohaboration with f« 5 fb^^ data at the 7 TeV run [74] and by the CMS 
collaboration with ~ 5 + 12 fb~^ data at the 7 TeV and 8 TeV runs [75]. Upper limits on 
the production times decay rates of these processes (which, unfortunately, have not given 
by the collaborations) have been set and they can be turned into constraints on the MSSM 
parameter space which, in the Higgs sector, corresponds to the [tan/3,M^] plane. 

In Fig. 7, we display the sensitivity of the CMS <I> — t- rr analysis in the [tan (3, Ma\ 
plane. The excluded region, obtained from the observed limit at the 95%CL is drawn 
in light blue. The solid line represents the median expected limit which turns out to be 
weaker than the observed limit. As can be seen, this constraint is extremely restrictive 
and for values Ma ^ 250 GeV, it excludes almost the entire intermediate and high tan/3 
regimes, tan/3 > 5. The constraint is of course less effective for a heavier pseudoscalar 
Higgs boson, but even for Ma = 400 GeV the high tan/3 > 10 region is excluded and one 
is even sensitive to large values Ma ~ 700 GeV for tan /3 > 50. 

There are, however, some caveats to this exclusion limit as discussed in section 3.1. The 
first one is that there is a theoretical uncertainty that affects the Higgs production cross 
section times decay branching ratios which is of the order of ±25% when the gg ^ ^ and 
66 —7- $ cross sections are combined. If this theoretical uncertainty is included when setting 
the limit in the [tan/3, M^] plane, as shown by the dashed contours around the expected 
limit in Fig. 7, the constraint will be slightly weaker as one then needs to consider the 
lower value of the (T(pp— >•$) xBR(<I> — )• t^t~) rate predicted by theory. 

The second caveat is that the CMS (and ATLAS) constraint has been given in a specific 
benchmark scenario, the maximal mixing scenario with the choice Xt/Ms = \/6 and the 
value of the SUSY scale set to Ms = 1 TeV; the other parameters such as the gaugino 
and higgsino masses and the first /second generation fermion parameters that have little 
impact can be chosen as in eq. (2.7). However, as was previously argued, the pp —>•$—)■ rr 
cross section times decay branching fraction is very robust and, hence, the exclusion limit 
is almost model independent. It is altered only very mildly by the radiative corrections in 
the MSSM Higgs sector, in particular by the choice of the parameters Ms and Xt (this is 
especially true if these parameters are to be traded against the measured values of M^). 

In fact, the exclusion limit in Fig. 7 can be obtained in any MSSM scenario with the 
only assumption being that SUSY particles are too heavy to affect a{pp — )■ $) x BR($tt) 
by more than 25%, which is the estimated theoretical uncertainty. Even in the case of light 
SUSY particles, it is very hard to make that stop/sbottom squarks contribute significantly 
to the gg ^ H production processes, or to have a significant A;, correction to the ^bb 
coupling which largely cancels out as indicated by eq. (2.12), or to have a substantial 
change of the ^ ^ tt fraction due to light SUSY particles that appear in the decays. 

Thus, the limit for the pp — t- t~^t~ searches is robust with respect to the SUSY param- 
eters and is valid in far more situations and scenarios than the "MSSM MF^^^ scenario" 
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that is usually quoted by the experimental collaborations. We thus suggest to remove this 
assumption on the benchmark scenario (in particular it adopts the choice Ms = 1 TeV 
which does not allow low tan/3 values and which starts to be challenged by direct SUSY 
searches), as the only relevant assumption, if any, should be that we do not consider cases 
in which the SUSY particles are too light to alter the Higgs production and decay rates. 
This is a very reasonable attitude if we are interested mainly in the Higgs sector. 
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Figure 7: The [tan/3, MA] plane in the MSSM in which the pp -> H/A -^ t+t^ (light blue) and 
t -^ bH^ -^ hTv (dark blue) search constraints using the CMS data are included (observed limits). 
The solid contour for the pp — >■ tt mode is for the median expected limit and the two dashed ones 
are when the QCD uncertainties on the rates are included. 

Another constraint on the MSSM Higgs sector^*^ is the one from charged Higgs searches 
in the H~ — t- tv final states with the H^ bosons produced in top quark decays, t — t- H^h — )• 
Tuh. Up to now, the ATLAS and CMS collaborations have released results only with the 
Ri 5 fb^^ collected at ^/s = 7 TeV [76,77]. We have also delineated in Fig. 7 the impact 
on the [tan /3, MA] parameter space of the CMS 95%CL observed limits in this channel. 

As can be observed, the constraint is effective only for values Ma ^ 150 GeV which 
correspond to a light H~^ state that could be produced in top quark decays. The search is 
sensitive to the very high tan /3 region which is completely excluded by the tt search, that 
is performed with much more data though. However, even if the comparison is made for 
the same amount of data, he pp — )• <1> — )• rr search is by far more sensitive. 

Note that contrary to the pp — > t'^t" case, the limits at high tan/3 from the process 
pp —7- tt with t — )• bH~^ —7- bTu might be more model dependent. Indeed, while SUSY decays 
might not be important as the small Mjj± value leaves little room for light sparticles 
(and the high tan/3 values would suppress these decays anyway), the effect of the A;, 
corrections might be larger as there is no cancellation between production and decay rates. 
Nevertheless, the H^ limit is effective only for Ma ^150 GeV values excluded by the h 
data. We keep this H^ constraint though, as it is also valid in two-Higgs doublet models. 



^'^A search has also been performed by the CMS coUaboration based on the 7 TeV data in the channel 
pp — >■ $66 — >■ bbbb [57] . This search is much less sensitive than the tt search even if one extrapolates the 
expected limits to the same amount of data. We will thus ignore it in our study. 



22 



4.3 Extrapolation to the low tan/5 region and the full 7+8 data 

A very important remark is that in our version of the constraints in the [tan/?, Ma] plane of 
Fig. 7, we have removed the region exchided by the bomid on the h mass, Mh ^114 GeV 
from negative Higgs searches at LEP2, which is also usuahy displayed by the experimental 
collaborations. In the usual benchmark scenario, this constraint excludes the entire low 
tan/3 regime, tan/3 < 3, and at low A-Ia ~ 100 GeV, tan/3 values up to tan/3 ~ 10. 

A first reason for removing the "LEP exclusion" region is that it is now superseded 
by the "observation" constraint 123 GeV < M^, < 129 GeV (once the theoretical and 
experimental uncertainties are included) which is by far stronger. In fact, as was discussed 
in Ref . [29] , if the benchmark scenario with Ms = 1 TeV and maximal stop mixing is to be 
adopted, the entire range tan/3 < 5 and tan/3 > 20 for any Ma value would be excluded 
simply by requiring that 123 GeV < Mh < 129 GeV (and the excluded regions would be 
completely different for other Ms and Xt values as also shown in Ref. [29]). 

A second reason is that the LEP2 Mh constraint and even the constraint Mh > 123 
GeV can be simply evaded for any value of tan/3 or Ma by assuming large enough Ms 
values as discussed in section 2.1. This will then open the very interesting low tan/3 region 
which can be probed in a a model independent way by Higgs search channels involving the 
H, A, H^ bosons, including the t — ?• hH^ — t- hrv channel discussed previously. 

Indeed, the branching fraction for the decay t — t- hH'^ is also significant at low tan/3 
values, when the component of the coupling gtbH+ that is proportional to rrii/tan/3 becomes 
dominant. On the other hand, the branching fraction for the decay H — > tv stays close 
to 100%. Hence, the rates for pp — )■ ti with t — )• hH^ — )• hru are comparable for tan/3 w 3 
and tan/3 ~ 30 and the processes can also probe the low tan/3 region. This is exemplified 
in Fig. 8 where the t — )■ hH^ CMS median expected and observed limits obtained with the 
7 TeV data are extrapolated to the low tan/S region. As can be seen, the region tan/3 < 2 
is excluded for M^i < 140 GeV (this region can also be probed in the -ff"*"— )-cs mode). 
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Figure 8: The [taii/3, MA] plane in the MSSM in which the pp -^ H/A -^ t+t^ (Uglit blue) and 

t — >■ 67J+ -^ hrv (dark blue) observed limits using the CMS data are extrapolated to low tan/3. 
The solid contours in the tt and tv cases are for the expected limits. 

In fact, as is shown in the lower part of Fig. 8, even the channel pp —t- $ —t- rr is useful 
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at low tan/3. Indeed, for tan/3 values close to unity, while the 66 — >• <I> process becomes 
irrelevant, the cross sections for the gg ^- ^ process becomes very large, the reason being 
that for tan/3 ~ 1 the couplings g^tt oc mf/tan/3 are significant and the dominant top 
quark loop contribution becomes less suppressed compared to the SM. On the other hand, 
at least in the case of the pseudoscalar A, the branching ratio for the t^t^ decay stays 
significant for M^ values up to the ti threshold as shown in Fig. 4. Hence, the production 
times decay rate for gg ^ A ^ tt stays large and the CMS search limit is eff'ective and 
excludes tan /3 values close to 1 , for pseudoscalar masses up to Ma ~ 350 GeV. 

One would get a better feeling of the power of these constraints at low tan/3 values 
(and in the charged Higgs case also at high tan/3), if the present limits in the pp — )■ tt and 
t —7- bH~^ — )• bru channels are extrapolated to the full set of data collected in the 2011 and 
2012 LHC runs. This is shown in Fig. 9 where the median expected CMS limits in the two 
search channels are extrapolated to an integrated luminosity of 25 fb~^, assuming that the 
limits simply scale like the square-root of the number of events. 

The gain in sensitivity is very significant in the H case as the gap between the present 
CMS limit with the ~ 5 fb~^ of the 7 TeV data and the expected limit with the additional 
~ 20 fb~^ data at 8 TeV is large (there is an additional increase of the pp — )■ ti production 
cross section from y/s = 7 TeV to 8 TeV). In the case of the pp — )■ tt channel, the increase 
of sensitivity is much more modest, not only because the gap from the 17 fb^^ data used 
in the latest CMS analysis and the full 25 fb~^ data collected up to now is not large but, 
also, because presently the observed limit is much stronger than the expected limit. 
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Figure 9: The [tan /J, MA] plane in the MSSM in which the pp -^ H/A - 
bri' CMS expected limits are extrapolated to the full 7+8 TeV data with 
observed limits are still shown in blue. 
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Hence, these interesting low tan/3 areas that were thought to be buried under the LEP2 
exclusion bound on M^ are now open territory for heavy MSSM Higgs hunting. This can 
be done not only in the two channels pp — t- t~^t~ and t — ?• bH^ — ?• bTi/ above (and which 
were anyway used at high tan/3) but also in a plethora of channels that have not been 
discussed before (or at least abandoned after the LEP2 results) and to which we turn now. 
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5. Heavy Higgs searches channels at low tan/3 

We come now to the main phenomenological issue of this paper: the probe at the LHC of 
the low tan/3 region for a not too heavy pseudoscalar A state^^. We stress again that this 
region can be resurrected simply by allowing a large SUSY scale Ms which removes the 
LEP2 Mh>lU GeV constraint (and now the LHC mass constraint Mh ^ 126 GeV). We 
show that several channels discussed in the case of a high mass SM Higgs or in scenarios 
beyond the SM can be used for the search of the MSSM H, A and H^ bosons. 

5.1 The main search channels for the neutral H/A states 

5.1.1 The H -^ WW, ZZ channels 

These are possible only for the heavier H boson (because of CP invariance there are no 
VV couplings for A) with masses below the ti threshold where the branching ratios for the 
decays H — t- WW and H — )• ZZ are significant; see Fig. 4. The H — >■ WW process is 
particularly useful in the region 160 < Adn ^180 GeV where the branching ratio is close 
to 100%. In both cases, the gg ^ H production process can be used but, eventually, vector 
boson fusion can also be relevant at the lowest tan/3 and Mh possible values. 

The search modes that are most useful at relatively low Mh values, should be the 
pp ^ H ^ ZZ — )■ 4^^ and pp ^ H ^ WW — ?• 2i2i' channels that have been used to 
observe the SM-like light h boson (as the mass resolution of the H — >• WW channel is rather 
poor, one has to subtract the observed signal events in the low mass range, Mh ^ 160 GeV) 
and to exclude a SM-like Higgs particle with a mass up to 800 GeV [69,71]. When the 
two processes are combined, the sensitivity is an order of magnitude larger than for the 
SM Higgs for masses below 400 GeV and one can thus afford a substantial reduction of the 
couplings gntt and gnw which should allow to probe tan /? values significantly higher than 
unity^^. At high H masses, Mh ^ 300 GeV, one could also add the pp ^ H ^ ZZ — >• 
2^2^^, 21^2^, 2^2z/ and pp ^ H ^ WW — )• lv2q channels to increase the statistics, as done 
in a recent study by the CMS collaboration [81]. 

There is one difference with the SM Higgs case though. While in the SM, the Higgs 
particle has a large total width at high masses as a result of the decays into longitudinal 
W/Z bosons which make it grow as Mfj , the MSSM H boson remain narrow as the 
coupling gnw is suppressed. In fact, all MSSM Higgs particles will have a total width 
that is smaller than ~ 3 GeV for tan /3 ~ 3 and masses below 500 GeV. The smaller total 
width in the MSSM can be rather helpful at relatively high H masses as, for instance, it 
allows to suppress the continuum ZZ background by selecting smaller bins for the invariant 



^^This issue has been discussed in the past and a summary can be found in section 3.3.2 of Ref. [7]. It 
has been also addressed recently in Ref. [53] (where, in particular, a feasibility study of the H ^t hh and 
A — >■ hZ modes at ^/a — 14 TeV is made) and in talks given the last months by one of the authors [79] . 
Recent analyses of heavier MSSM Higgsses at intermediate and high tan /3 can be found in Ref. [80] . 

^^The ATLAS collaboration has recently analyzed heavy H production in a two-Higgs doublet model in 
the channel H — >■ WW — >■ eufiv with 13 fb~^ data collected at yG = 8 TeV [82]. Unfortunately, this analysis 
cannot readily be used as the limit on the cross section times branching fraction has not been explicitly 
given and the results are displayed in terms of cos(a) (and not cos(/3 — a) which would have corresponded 
to the HWW coupling) which does not allow an easy interpretation in the MSSM. 
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mass of the ZZ system in the signal events. Issues like the interference of the signal and 
the gg — )■ W backgrounds will also be less important than in the SM. 

5.1.2 The H/A -^ tt channels 

This search channel has not been considered in the case of the SM Higgs boson for two 
reasons [4]. The first one is that for M^^^ > 350 GeV, the i?sM — ^ WW,ZZ channels 
are still relevant and largely dominate over the H^m -^ tt decay channel which has a 
branching fraction that is less than 20% in the entire Higgs mass range (the reason being 
again that the partial widths for -ffsM — ^ VV grow as Mfj while for Hsm — ^ ti its grows 
only like Mh^m ) • '^^^ other reason is that the continuum tt background was thought to be 
overwhelmingly large as it had to be evaluated in a large mass window because of the large 
Higgs total width (in addition, the events from -ffsM — ^ tt produce a dip-peak structure in 
the gg — )• ti invariant mass spectrum that was unobservable for a large total width) . 

The situation in the MSSM is very different. First, as mentioned previously, the total 
width for heavy H and A states are much smaller, less than < 20 GeV for any tan/3 > 1 
value for Mh^a ^ 500 GeV and grow (almost) linearly with the Higgs masses beyond this 
value. One can thus integrate the ti continuum background in a smaller invariant mass 
bin and significantly enhance the signal to background ratio. A second feature is that 
contrary to the SM case, the branching ratios for the H/A — t- ti decays are almost 100% 
for tan/3 < 3 as soon as the channels are kinematically open (this is particularly true for 
A where even below the threshold, the three-body decay A — )■ tt* — t- tbW is important). 

The only disadvantage compared to the SM is that the production cross section could 
be smaller. In the MSSM, the only relevant process in the low tan/3 regime for M<j, > 350 
GeV \s gg ^ ^ with the dominant (almost only) contribution being due to the top quark 
loop. The latter is suppressed by the square of the coupling g^tt oc l/tan/3 if tan^S is not 
close to unity. However, in the MSSM, one has to add the cross sections of both the H 
and A states. In addition, the loop form factors in the pseudoscalar A and scalar H/Hsm 
cases are different and, as can been seen from Fig. 3, the gg ^ ^ cross section is larger in 
the pseudoscalar Higgs case when the same top Yukawa coupling is assumed. 

In toto, the situation for H/A — )• ti will certainly be more favorable for the MSSM at 
low tan /3 than in the SM. While there was no search for the SM Higgs in this channel, the 
ATLAS [83] and CMS [84] collaborations have looked for heavy resonances (such as new 
Z' gauge bosons in extended gauge models or Kaluza-Klein excitations in scenarios with 
extra space-time dimensions) that decay into ti pairs with the data collected at the 7 TeV 
run. The lepton+jets final state has been studied in the topology where the top quarks 
are highly boosted which allows a good discrimination from the continuum ti background 
[85] (the ATLAS and CMS collaborations searches assume resonance masses Mu ^ 700 
GeV to benefit from this topology). Limits on the cross sections times branching ratios 
have been set, corresponding to roughly au ~ 0.7 pb for a resonance with a mass of 1 TeV 
and a narrow width, Ttt ~ lO^^M^j (which is more or less the case of the MSSM H/A 
states at tan/3 ~ 3). A lower (higher) cross section is needed at larger (smaller) resonance 
mass when the top quarks are (not) sufficiently boosted and, at Mft ~ 500 GeV, one needs 
ctt ~ 3 pb which approximately corresponds to an increase with 1/M^^. 
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5.1.3 The A ^ Zh channel 

As discussed earlier, the gg ^ A production cross section is very large at low tan /? values: 
it is higher than for the SM Higgs boson at tan/3 = 1 (as the form factor for the ggA 
amplitude is larger than in the scalar Higgs case) and is suppressed only by a factor 
g\^^ oc l/tan^/3. On the other hand, in the range M^ + Mz ^ Ma ^ 2mt, the branching 
ratio for the decay A — )■ hZ is large for tan/3 ~ 3 and largely dominant for tan/3 ~ 1. 
In the mass window Ma = 210-350 GeV, the production times decay rate for the process 
gg ^ A ^ hZ should be thus very high in the low tan/3 region. 

The hZ final state has been searched for in the SM in the Higgs-strahlung process, 
qq ^ Z* -^ Zh with the Z boson decaying into leptons or neutrinos, Z — t- i!.'^i.~,vi> and 
the h boson decaying into 56 final states [72, 73]. The significance of the signal is strongly 
increased by looking at boosted jets when the Higgs has a large transverse momentum 
[86]. In the CMS analysis with 17 fb^^ of the 2011 and 2012 data [73], a signal strength 
/ifefe « 1.5 has been found in the Z ^ uv and Z — )■ I'^i" channels with a large error bar. 
Very roughly, one can assume that the additional events from the A — > Zh channel should 
be observed if they exceed this sensitivity when extrapolated to include the full 2012 data. 

One should note that the information from the pp — )■ Zh search in the SM provides 
only a lower limit for the sensitivity as in the present case one can benefit from the fact 
that the invariant mass of the four fermion final state (without neutrinos) which should 
peak at the value Ma will further suppress the continuum background, in particular the 
Z + bb events. However, as h is originating from the decay of the state A which should 
not be very heavy, it has not enough transverse momentum to strengthen the boosted jet 
techniques that allow to isolate the h ^ bb signal from the QCD background. 

5.1.4 The H — > hh channel 

The channel pp ^ H ^ hh is similar to A — )■ hZ: it has very large production rates in the 
low tan/3 regime in the mass range 250 GeV < Mh ^ 350 GeV when the decay channels 
i7 — 7- /i/i is kinematically open and the H ^ ti mode is closed; the gg ^ H cross section 
should be substantial in this area of the parameter space. 

If the dominant h ^ bb decay is considered, the signal topology has some similarities 
with that of the process gg ^ bb + A/H which was discussed here as being one of the 
main MSSM Higgs processes at high tan/3 and searched for by the CMS collaboration 
with the 7 TeV data [57]. However, the kinematical behavior is very different and in the 
signal events, one can use further constraints, M^fe ~ M^ and Mfe^^f, ~ Mh (see Ref. [53] 
where a characterization of this channel has been made). In fact, the H ^ hh channel has 
more similarities with double production of the SM-like Higgs boson, gg — )■ /i/i, which is 
considered for the measurement of the Higgs self-coupling a the 14 TeV LHC with a high 
luminosity. This process has been revisited recently [87] and it has been shown that the 
final state channels bbrr and 6677 would be viable at y/s = 14 TeV and C > 300 fb^^. 
Because the /i — t- 77 decay is too rare, only the first process could be considered at -^s = 8 
TeV with 25 fb~^ data. Note that here again, one could use the reconstructed H mass 
constraint, Mh = M^^, to further suppress the continuum background. 
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5.2 Expectations for the LHC with 25 fb ^ data 

It is obvious that a truly rehable estimate of the sensitivity on the heavy neutral MSSM 
Higgs bosons in the various channels discussed before can only come from the ATLAS and 
CMS collaborations. We will nevertheless attempt in this section to provide a very rough 
estimate of the achievable sensitivities using present searches conducted for a heavy SM 
Higgs and in beyond the SM scenarios. The very interesting results that could be obtained 
would hopefully convince the experimental collaborations to conduct analyses in this area. 
Following the previous discussions, our working assumptions to derive the possible 
sensitivities in the various considered search channels are as follows: 

- iJ — )■ WW, ZZ: we wiU use the recently pubhshed CMS analysis of Ref. [81] that 
has been performed with the ~ 10 fb~^ data collected in the 7+8 TeV runs and in which 
all possible channels H ^ ZZ ^ At, 2i2u, 2£jj, 2ujj and H -> WW -^ 2i2u, (.vjj have 
been included and combined. In the entire range Mh = 160-350 GeV, where the SM Higgs 
boson almost exclusively decays into WW or ZZ states, we will assume the cross section 
times decay branching ratio upper limit that has been given in this CMS study, 

~ H/A -f ti: we make use of the ATLAS [83] and CMS [84] searches at ^/s = 7 TeV for 
new Z' or Kaluza-Klein gauge bosons that decay into ti pairs in the lepton+jets final state 
topology. Considering a small total width for the resonance, limits on the cross sections 
times branching ratio of ~ 6, 3 and 0.75 pb for a resonance mass of, respectively, 350, 500 
and 1000 GeV are assumed. This is equivalent to a sensitivity that varies with l/M^f that 
we will optimistically assume to also cover the low mass resonance range. 

- A —7- hZ: we will use the sensitivity given by ATLAS [72] and CMS [73] in their 
search for the SM Higgs-like strahlung process pp — )• hZ with h ^ bb and Z — )■ ii, vV, 
cr/cjS^ = 2.8 with 17 fb~^ data at y/s = 7 + 8 TeV (we have included the error bar). 
This should be sufficient as, in addition, we would have on top the constraint from the 
reconstructed mass in the Ubb channel which is not used in our analysis. 

- H ^ hh: we use the analysis of the process gg — )■ hh in the SM performed in Ref. [87] 
for the 14 TeV LHC that we also scale down to the current energy and luminosity. The 
final state bbrr final state will be considered, with the assumption that the cross section 
times branching ratio should be larger than a x BR ~ 500 fb for illustration. 

The results are shown in Fig. 10 with an extrapolation to the full 25 fb~^ data of the 
7+8 TeV LHC run. Again, we assumed that the sensitivity scales simply as the square root 
of the number of events. The sensitivities from the usual H/A — )• t'^t" and t — )• bH'^ — )• bri' 
channels are also shown. The green and red areas correspond to the domains where the 
H — >■ VV and H/A — )• ti channels become constraining with the assumptions above. The 
sensitivities in the H ^ hh and A — )■ hZ modes are given by, respectively, the yellow and 
brown areas that peak in the mass range MA = 250-350 GeV visible at very low tan/3 
values. We refrain from extrapolating to the LHC with 14 TeV cm. energy. 

The outcome is impressive. These channels, in particular the H — t- VV and H/A — t- ti 
processes, are very constraining as they cover the entire low tan/3 area that was previously 
excluded by the LEP2 bound up to M^ ~ 500 GeV. Even A — )• hZ and H ^ hh are visible 
in small portions of the parameter space at the upgraded LHC. 



28 



a 




200 



400 

Ma [GeV] 



600 



800 1000 



Figure 10: The estimated sensitivities in the various search channels for the heavier MSSM Higgs 
bosons in the [tan/3, Myi] plane: H/A -^ t+t^ (light blue), t -^ H+b -^ rvb (dark blue), H — > 
WW + ZZ (green), H/A -^ tt (red), A ^ hZ (brown) and H ^ hh (yellow). The projection 
is made for the LHC with 7+8 TeV and the full 25 fb~^ of data collected so far. The radiative 
corrections are such that the lightest h mass is M^ = 126 GeV. 

5.3 Remarks on the charged Higgs boson 

We close this discussions with a few remarks on the charged Higgs boson case. First of all, 
the production rates are very large only for Mjj± < 170 GeV when the H^ state can be 
produced in top decays. In this case, the decay channel H — t- tv is always substantial and 
leads to the constraints that have been discussed earlier and which are less effective than 
those coming from H/A — t- tt searches at high tan/3. In the low tan/3 region, two other 
channels can be considered: H^ — t- cs that has been studied by the ATLAS collaboration 
in a two-Higgs doublet model with the 7 TeV data [88] and H^ — t- cb. The branching ratio 
for the latter channel is significant for tan/3 < 3 and has been obtained by assuming the 
same CKM angles as in the SM, in particular Vcb ~ 0.04 [34]. This channel, if observed 
would thus allow to check some of the CKM matrix elements in the charged Higgs sector. 

Finally, the processes t — >• H~^b at low mass and pp — )• btH^ at high mass with 
H^ —7- Wh can have large rates at sufficiently low tan/3. The cross section times branching 
fraction is displayed in Fig. 11 in the [tan/3,Af^] plane for a 14 TeV cm. energy. Shown 
are the contours with a x BR = 1.10 and 30 fb which, for a luminosity of 300 fb~^ would 
correspond to a small number of events. We will not perform an analysis for this particular 
final state. We simply note that the final state topology, pp — )• tbH — )• tbWh resembles 
that of the pp — )■ tth process that is considered as a means to measure the htt Yukawa 
coupling and which is considered to be viable at 14 TeV with a high luminosity. 

Hence, even for the charged Higgs bosons, there are interesting search channels which 
can be considered if the low tan /3 region is reopened. 
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Figure 11: The production cross sections times decay branching ratio at the LHC with ^/s = 14 
TeV for the process pp — ^ tbH^ + ibH^ with H^ — >■ hW^ in the [tan/?,M^] plane. The contours 
are for some (rough estimates of) limiting values of cr x BR. 

6. Conclusions 

After the observation of the 126 GeV SM-hke Higgs boson by the ATLAS and CMS 
collaborations, the next challenge at the LHC should be to search for new phenomena 
beyond the SM. This can be done not only by refining the precision determination of the 
properties of the observed Higgs particle to pin down small deviations of its couplings from 
the SM expectations, but also by looking for the direct production of new states. 

In this paper, we have considered the production of the heavier H, A and H bosons of 
the MSSM at the LHC, focusing on the low tan/3 regime, tan/3 < 3-5. We have first shown 
that this area of the MSSM parameter space, which was long thought to be excluded, is 
still viable provided that the SUSY scale is assumed to be very high, Ms ^ 10 TeV. For 
such Ms values, the usual tools that allow to determine the masses and couplings of the 
Higgs and SUSY particles in the MSSM, including the higher order corrections, become 
inadequate. We have used a simple but not too inaccurate approximation to describe the 
radiative corrections to the Higgs sector, in which the unknown scale Ms and stop mixing 
parameter Xt are traded against the measured h boson mass, M^ ~ 126 GeV. One would 
then have, to a good approximation, only two basic input parameters in the MSSM Higgs 
sector even at higher orders: tan/3 and M^, which can take small values, tan/3 w 1 and 
Ma = 0(200) GeV, provided that Ms is chosen to be sufficiently large. 

In the low tan /3 region, there is a plethora of new search channels for the heavy MSSM 
Higgs bosons that can be analyzed at the LHC. The neutral H/A states can be still be 
produced in the gluon fusion mechanism with large rates, and they will decay into a variety 
of interesting final states such as H ^ WW, ZZ, H — ;• hh, H/A -^ ti, A ^ hZ. Interesting 
decays can also occur in the case of the charged Higgs bosons, e.g. H^ — t- hW, cs, cb. These 
modes come in addition to the two channels H/A — t- t^t^ and t — >• bH'^ — t- bri/ which are 
currently being studied by ATLAS and CMS and which are very powerful in constraining 
the parameter space at high tan/3 values and, as is shown here, also at low tan/3 values. 
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We have shown that already with the current LHC data at ^/s = 7+8 TeV, the area 
with small tan /3 and Ma values can be probed by simply extrapolating to the MSSM Higgs 
sector the available analyses in the search of the SM Higgs boson at high masses in the WW 
and ZZ channels and the limits obtained in the ti channels in the search for high-mass new 
gauge bosons from extended gauge or extra-dimensional theories. The sensitivity in these 
channels will be significantly enhanced at the 14 TeV LHC run once 300 fb~^ data will be 
collected. In the absence of any signal at this energy, the [tan/3, Ma\ plane can be entirely 
closed for any tan /? value and a pseudoscalar mass below Ma ~ 500 GeV. Additional and 
complementary searches can also be done in the charged Higgs case in channels that have 
not been studied so far such as H^ — t- Wh but we did not analyze this issue in detail. 

Hence, all channels that have been considered for the SM Higgs boson in the high mass 
range, plus some processes that have been considered for other new physics searches, can 
be recycled for the search of the heavier MSSM Higgs bosons in the low tan /3 regime. For 
instance, many of these MSSM Higgs processes could benefit from the current searches of 
multi-lepton events with missing energy in SUSY theories. As in all channels we have W, Z 
and additional h bosons in the final states, multileptons and missing energy are present in 
most of the topologies. One could then use the direct searches for SUSY particles such as 
charginos and neutralinos to probe also the MSSM heavier Higgs states. 

All this promises a very nice and exciting program for Higgs searches at the LHC in 
both the present and future runs. One could then cover the entire MSSM parameter space: 
from above (at high tan/3) by improving the H/A — t- tt searches but also from below (at 
low tan/3) by using the WW^ ZZ^ tt, .. searches. The coverage of the [tan/3. Ma] plane will 
be done in a model independent way, with no assumption on Ms and possibly on any other 
SUSY parameter ^^. The indirect information from the lighter Higgs mass will be included 
as well as the information from the Higgs couplings, as the sensitivity regions cover also 
that which are excluded from the measurement of the h properties at the LHC. 

One can of course use these channels in other extensions of the SM. An example would 
be SUSY extensions beyond the MSSM where M^ can be made large enough without 
having large Ms values; this is the case of the NMSSM where the maximal Mh value can 
be obtained at tan/3 ss 2 [89]. Another example would be a non-SUSY two-Higgs doublet 
model where there is more freedom in the parameters space and all channels analyzed here 
and even some more could be relevant; discussions along these lines have already started 
[82,90]). The numerous search channels discussed in this paper might allow to probe in a 
more comprehensive manner the extended parameter space of these models. 

Acknowledgements: AD thanks the CERN theory division for the kind hospitality of- 
fered to him. This work is supported by the ERC Advanced Grant Higgs@LHC. 



^^This approach is orthogonal to that of Ref. [32] in which specific benchmark scenarios with fixed SUSY 
parameters (which might need to be updated soon) are proposed. We note that for all the proposed 
benchmarks scenarios [32], the SUSY scale is fixed to Ms = 1 or f.5 TeV which excludes the low (and 
possibly intermediate) tan /3 regime and, hence, the possibility of discussing the processes analysed here. 
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